Isolated resonances in conductance fluctuations and hierarchical states 
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We study the isolated resonances occurring in conductance fluctuations of quantum systems with 
a classically mixed phase space. We demonstrate that the isolated resonances and their scattering 
states can be associated to eigenstates of the closed system. They can all be categorized as hier- 
archical or regular, depending on where the corresponding eigenstates live in the classical phase 
space. 



I. INTRODUCTION 

The classical dynamics of a scattering system is re- 
flected in the transport properties of its quantum me- 
chanical analog. A prominent example in quantum chaos 
are the universal conductance fluctuations exhibited by a 
scattering system with classically completely chaotic dy- 
namics [y. Generic systems, however, are neither com- 
pletely chaotic nor intcgrable, but show chaotic as well 
as regular motion The chaotic dynamics is strongly 
influenced by the presence of islands of regular motion, 
in particular, one finds a trapping of chaotic trajectories 
close to regular regions with trapping times distributed 
according to power laws . The semiclassical analysis re- 
vealed that conductance fluctuations of generic scattering 
systems have corresponding power-law correlations 
and most interestingly that the graph conductance vs. 
control parameter is a fractal 0. Fractal conductance 
fluctuations have indeed been observed experimentally 
in semiconductor nanostructures ^,0] as well as numeri- 
cally §. 

Surprisingly, for the cosine billiard j^jltj, a system 
with a mixed phase space and power-law distributed 
classical trapping times, a recent numerical study did 
not show fractal conductance fluctuations (If). Instead, 
sharp isolated resonances were found with a width distri- 
bution covering several orders of magnitude. Only about 
one third of them can be related to quantum tunneling 
into the islands of regular motion |jl2), while the rest 
remained unexplained. It was later shown that conduc- 
tance fluctuations for mixed systems should in general 
show fractal fluctuations on large scales and isolated res- 
onances on smaller scales |l3j . The isolated resonances in 
the scattering system were conjectured to be related to a 
subset of eigenstates of the closed system, namely hier- 
archical states jl4j living in the chaotic component close 
to the regular regions and regular states living within the 
islands of regular motion Jl^j . This type of behavior was 



obtained for a quantum graph which modeled relevant 
features of a mixed phase space |0 . 

The purpose of the present paper is to establish the ori- 
gin of all isolated resonances for a system with a mixed 
phase space. To this end, we study the cosine billiard for 
suitable parameters in a three-fold way: (i) as a quan- 
tum scattering system, (ii) as a closed quantum system, 
and (iii) its classical phase-space structures. We find that 
the resonances have scattering states and corresponding 
eigenstates of the closed system that live in the hierar- 
chical and regular part of phase space. The number of 
resonances of each type is directly related to the cor- 
responding volumes in the classical phase space. Each 
resonance width is quite well described by the strength 
of the corresponding eigcnfunction at the billiard bound- 
ary. Exceptions are shown to arise from the presence of 
avoided crossings in the closed system. It is demonstrated 
that the simultaneous appearance of fractal conductance 
fluctuations and isolated resonances, as observed in a 
quantum graph model p3| , would for our system with a 
mixed phase space require much higher energies. These 
are currently computationally inaccessible. 

In the following section, we define the model we use to 
study the relation between the scattering resonances and 
the eigenstates of the corresponding closed system. Our 
main results on the classification of resonant scattering 
states and corresponding eigenstates of the closed system 
into hierarchical and regular are presented in Sees. Ill 
and |y|. The role of partial transport barriers is analyzed 
m Sec. [V|. In Sec. |VT| we discuss the effect of avoided 



crossings on the assignment of resonances of the o pen 
billiard to eigenstates of the closed system and Sec. VII 
gives a summary of the results. Finally, the Appendix 
contains some details of the numerical methods employed 
in the present work. 
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II. THE MODEL 

We study the cosine billiard J|Jl(]] , either closed or with 
semi-infinite leads attached. The boundaries of the bil- 
liard are hard walls (i.e. Dirichlet boundary conditions) 
at y — and 



/ \ ttt M ( , [2-kx\ 



(1) 



for < x < L (see Fig. |l]a) . In the open billiard two semi- 
infinite leads of width W are attached at the openings at 
x = and x = L, while in the closed billiard the open- 
ings are closed by hard walls. The classical phase-space 
structure can be changed by varying the ratios W / L and 
AI/L. For M/L = the dynamics is integrable and for 
example for M/L =1/2 and W/L = 1 the dynamics ap- 
pears to be ergodic (at least the islands of regular motion, 
if any, are very small) || . 

In the present work, we use the same parameters as 
in Ref. [jll|, namely W/L = 0.18 and M/L = 0.11, for 
which the I- and M-shaped orbits depicted in Fig. |l|a 
are stable. The corresponding Poincare section is shown 
in Fig. [j]b. We use Poincare-Birkhoff coordinates (s,p), 
where s is the arclength along the upper boundary of the 
billiard with length L' « 1.029 L and p is the projection 
of the unit momentum vector after the reflection on the 
tangent in the point s. 

Quantum mechanically, for a given wave number Uf 
the number N of transmitting modes in a lead of width 
W is the largest integer with TV < kpW/ir. We measure 
energies in units of the energy Eq = h 2 n 2 / (2m,W 2 ) of the 
lowest mode in such a lead, i.e., E = h 2 'k 2 F / '(2mEo) = 
(UfW/tt) 2 > N 2 . The larger the number N of modes 
is, the more details of the classical phase space can be 
resolved by quantum mechanics. At the same time the 
computational effort increases with iV 4 and we compro- 
mise, as in Ref. |llj , on the case of N — 45 transmitting 
modes in the energy range E S [2026, 2100]. 



III. RESONANCES AND SCATTERING STATES 

Resonances in the scattering system, which have been 
observed as isolated features in conductance fluctua- 
tions |ll[ , were identified by isolated peaks in the Wigner- 
Smith time delay r of the system. The time delay is given 
by 



r = w Tr (SUS/dE), 



(2) 



where 2N is the dimension of the S-matrix. The calcu- 
lation of S and t was already outlined in Ref. [[II) and is 
presented in greater detail in Appendix |X| 

In Fig. we show the Wigner-Smith time delay r [in 
units of h/E = 2mW 2 /(hir 2 )} for E e [2026,2100]. The 
isolated resonances found in Ref. [^l) are clearly seen. 
Each isolated resonance E ICSt i has a Breit-Wigner shape 
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FIG. 1. (a) The cosine billiard with semi-infinite leads 
(short dashed lines) and hard walls for closing the system 
(dotted lines) for W/L = 0.18 and M/L = 0.11. Also 
shown are the two most prominent stable periodic orbits 
for these parameters (long dashed lines), (b) Poincare sec- 
tion of some regular and one chaotic orbit for the above pa- 
rameters in Poincare-Birkhoff coordinates p vs arclength s 
along the upper boundary of the billiard. A major island 
at (s,p) = (L'/2,0) around the elliptic I-shaped orbit and 4 
smaller islands surrounding the M-shaped orbit can be seen. 
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with TiTi = 2/N. Note that the heights and the cor- 
responding widths Ti of the individual resonances cover 
several orders of magnitude. 

In order to elucidate the nature of the resonances, we 
calculated the scattering states inside the open billiard. 
For a given configuration of waves incoming in both leads, 
the knowledge of the S'-matrix allows the determination 
of the outgoing waves and hence the wavefunction ampli- 
tudes at the openings of the billiard. Since the S'-matrix 
is defined between asymptotic, propagating modes, this 
procedure neglects the contribution of evanescent modes 
in the leads in the vicinity of the billiard. The wave- 
function amplitudes at the openings can then be used as 
boundary conditions for the solution of the Schrodinger 
equation inside the billiard. For the examples of scat- 
tering states presented below, we occupied the 10 top- 
most modes incoming from the left with equal ampli- 
tudes. Similar pictures were obtained for other boundary 
conditions. 

For the comparison of the scattering states with 
the classical phase-space structures we have calculated 
Husimi projections H sc (s,p). Similar to the case of 
closed billiards, see sec. 0, we define these by the pro- 
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FIG. 2. Wigner-Smith delay time r vs energy E. For each resonance a corresponding eigenstate of the closed system was 
found and the labels indicate whether it lives in the regular (r) or hierarchical (h) region of phase space. 



jection of the scattering state onto a coherent state on 
the upper boundary of the billiard, 



H sc (s, P ) = |(9 n ^|0™ hcrcnt (s, P ))| 2 



L' 



ds'd n ip*{s')e ikp{s '^ s) -^ k{s '- s)2 



(4) 



(5) 



with k = y/En/W. Here d n ip(s) = n(s) • VV>(q(s)) is the 
normal derivative of the scattering state on the upper 
boundary. n(s) is the normal vector and q(s) is the posi- 
tion of the boundary as a function of arc length s. Note 
that these Husimi projections are not normalized and are 
influenced by the openings over a range of a few Fermi 
wavelengths. Also they do not include the full billiard 
boundary and therefore no periodization of the coherent 
state has been used. 

As a first example, we present in Fig. || on the left the 
scattering state at an energy of approximately 2029.172, 
the center energy of the sharpest observed resonance. 
Obviously, the scattering state is associated with the I- 
shaped periodic orbit. The wavefunction amplitude is 
concentrated near the orbit and the Husimi projection 
lives predominantly inside the central stable island of 
the classical phase space. For comparison, we present in 
Fig. H on the right the scattering state at energy 2041.109. 
The width of the resonance at this energy is about 6 ■ 10 7 
times larger than the width of the sharpest resonance. 
Evidently, this resonance is not related to the stable is- 
lands in phase space. In contrast, by comparing with 
the superimposed KAM-tori of the Poincare section and 



a partial transport barrier surrounding the island hierar- 
chy (see Section [v]), we see that the Husimi projection 
lives in the hierarchical region between the islands and a 
partial transport barrier. 

As scattering states allow a great variability in the 
boundary conditions, e.g., the incoming modes, we do 
not use them for a detailed analysis of the isolated res- 
onances. Instead we consider the corresponding eigen- 
states of the closed system in the next section. 



IV. RESONANCES AND CORRESPONDING 
EIGENSTATES 

In this section we want to demonstrate that the iso- 
lated resonances of the conductance fluctuations and 
their scattering states have corresponding eigenstates of 
the closed billiard. In particular, we will show that all 
these eigenstates live in the hierarchical and regular part 
of phase space, as was conjectured in Ref. O]. This 
allows a labeling of all isolated resonances appearing in 
Fig. | 

For the closed system the eigenvalues and eigenfunc- 
tions are computed using the boundary element method, 
see, e.g. []15| and references therein. Since the cosine bil- 
liard is symmetric with respect to the axis x = L/2, 
the eigenstates have definite parity P = +, — . The ac- 
tual calculations are performed for the desymmetrized 
billiard with either Dirichlet or Neumann boundary con- 
ditions on the symmetry axis yielding the antisymmet- 
ric (P = — ) and symmetric (P = +) states, respec- 
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FIG. 3. Resonant scattering states (top row) and Husimi projections (bottom row) onto the classical Poincare section 
with KAM-tori (solid lines) and a partial transport barrier (dashed line). Two example at resonance energies 2029.172 (left) 
and 2041.109 (right) are shown. They live in the regular and the hierarchical region of phase space, respectively. For the 
representation of the scattering states a superposition with equal weight of the 10 topmost modes incoming from the left is 
shown. 



tively. We label the n-th eigenstate of parity P by n p . 
The mean level spacing A is determined by the area 
A = L(W + M/2) of the billiard using Weyl's formula 
A/E = (4:TTh 2 /2mA)/Eo = 0.176. 

We present in Fig. [| the two eigenstates correspond- 
ing to the scattering states shown in Fig. |^. For each 
state, we show the eigenfunction density |?/> n (q)| 2 and 
the corresponding Husimi representation H n (s,p) (see 
e.g., |l6|Jl7|| ). The state 5686~, displayed on the left of 
Fig. W differs in energy by about 0.01A from the sharpest 
observed resonance with energy 2029.172. Note that this 
energy difference is of the order of the accuracy to which 
our resonance energies and eigenvalues are calculated. 
On the right hand side of Fig. a hierarchical state is 
displayed. Its energy differs from the resonance at energy 
2041.109 by about 0.1A. This shift of the resonance en- 
ergy from the eigenenergy of the closed system is due to 
the opening of the system by attaching the leads. As for 
the scattering states, we have superimposed some KAM- 
tori onto the Husimi representations of Fig. ||. In addi- 
tion, a partial transport barrier surrounding the island 
hierarchy is shown (see next section). 

Now we want to associate all resonances of the scatter- 
ing system with width T at energy E res with an eigenstate 
of the closed billiard with energy E cv . We use a Husimi 
representation H n (s,p) on the Poincare section to deter- 
mine the region on which an eigenstate localizes. We 
introduce the quantity 



Vn = ds dpH n (s,p), (6) 
J-W J-l 

which integrates the Husimi distribution over the left 
boundary of the billiard (not shown in Fig. |J) with the 
normalization of the Husimi distribution chosen such that 
I-w 2 ds J^dp H n (s,p) = 1. This quantity gives an es- 
timate on how strongly a state of the closed system will 
couple to the leads in the scattering system and should 
be roughly proportional to T. This allowed us to find, 
for each of the 54 resonances with T < A/2, a state with 
E cv ss E Ies and with r/ w r. Figures ^ and ^ show the 
difference E cv — E ICS in units of the mean level spacing A 
and the approximate proportionality of rj and T, respec- 
tively. Clearly the larger differences appear for bigger T, 
but still a clear identification is possible (see section [y|) . 
This assignment also works the other way around, as of 
the 46 eigenstates with the smallest values of rj, we can 
identify 40 with isolated resonances, missing only the 6 
regular eigenstates quantized most deeply in the central 
island of phase space, as discussed below. 

For the 54 resonances with width T less than half a 
mean level spacing, we analyze the structure of the cor- 
responding eigenstates. We find that 17 states can be 
categorized as regular states, as their Husimi represen- 
tations live inside the five major stable islands in phase 
space. Of these states, 7 are associated with the I-shaped 
orbit and 10 with the M-shaped orbit. While we observe 
all states in the energy interval associated with the M- 
shaped orbit, 6 further eigenstates living near the center 
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of the central stability island are not resolved as reso- 
nances. As these are the innermost states in the island, 
we expect them to couple weaker to the leads and their 
resonance widths to be much smaller than the sharpest 
observed resonance. Apparently, these resonances are so 
narrow that we were not able to find them, given our cur- 
rent numerical accuracy, even knowing their approximate 
energy from the eigenvalues. 

The remaining 37 resonances are not related to regu- 



lar states, but the Husimi representations of their cor- 
responding eigenstates have large intensity in the region 
between the regular islands and the partial transport bar- 
rier and a much weaker intensity in the rest of the chaotic 
region. It should be noted that in the studied energy 
range accessible to our methods the wave length is of 
the order of the distance between regular islands and the 
partial transport barrier. Therefore the eigenstates ei- 
ther look like regular states living outside the island |I§|] 
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FIG. 5. Difference of eigenstate energy E cv and resonance 
energy E res in units of the mean level spacing A vs T/A. The 
deviations increase with F. 



FIG. 6. The strength 77 of an eigenstate at the left bound- 
ary vs the resonance width F of the corresponding resonance. 
An approximate proportionality can be seen. 
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or similar to scarred states on a hyperbolic orbit close 
to the island |is[| . For higher energies they would show 
the true properties of hierarchical states, i.e., being sim- 
ilar to a chaotic state, but restricted to the hierarchical 
region H]. We therefore classify these states as hierar- 
chical states. 

In Fig. H we have labeled the resonances by r and 
h according to our classification of the corresponding 
eigenstates as regular and hierarchical, respectively. This 
demonstrates that the origin of all isolated resonances are 
hierarchical or regular eigenstates of the closed system. 



V. PARTIAL TRANSPORT BARRIERS 

Classical transport in the chaotic part of phase space is 
dominated by partial barriers |20|-|24|]. They are formed 
by cantori as well as by stable and unstable manifolds. 
Such a partial transport barrier coincides with its it- 
erate, with the exception of so-called turnstiles where 
phase-space volume is exchanged between both sides of 
the partial barrier. We have constructed partial barriers 
using the methods described in Ref. (20). The fluxes <& 
are determined from the length I of the maximizing and 
minimax orbits, according to 



$ Tlkp | ^maximizing ^minimax | 



maximizing ''minimax 



\/w. 



(7) 
(8) 



Quantum mechanically, partial transport barriers with 
fluxes up to the order of h divide the chaotic part of phase 
space into distinct regions with chaotic and hierarchical 
eigenf unctions living mainly on one or the other side jl4j . 
We found that the partial barrier with smallest flux that 
surrounds the main island and the 4 neighboring islands 
can be constructed from the stable and unstable mani- 
folds of the period 4 hyperbolic fixed points. Each of its 
two turnstiles has for our largest energy E — 2100 a flux 
1.06?i. Further outside are many other partial barriers 
with only slightly bigger fluxes. As an example, we show 
in Figs, |3j and |4] the partial barrier constructed from an 
unstable periodic orbit with winding number 5/23, which 
is an approximant of the most noble irrational between 
winding numbers 1/4 and 1/5. It has a flux 1.65?i. 

A check on the validity of our identification of regular 
and hierarchical states is provided by a comparison of 
their numbers to the corresponding relative volumes in 
phase space. To this end, we calculate the volume of the 
tori associated with stable periodic orbits, V r , and the 
chaotic phase-space volume inside the partial transport 
barrier, Vh- We find V r and Vh to cover 5.9% and 8.5 % 
of the energy shell, respectively. From the total number 
of eigenstates in the energy interval, N — 426, we get 
for 23 (17+6) regular and 37 hierarchical states relative 
fractions of 5.4% and 8.7%, respectively, in good agree- 
ment with the volumes of the associated regions in phase 
space. 



The absence of fractal conductance fluctuations in 
this system has now a clear explanation: According to 
Ref. pj|] for fractal fluctuations to occur a hierarchy of 
partial transport barriers with fluxes larger than h must 
exist. For the studied energies we find that even the 
outermost partial barriers surrounding the hierarchical 
phase-space structure have fluxes of the order of h. This 
causes a quantum dynamical decoupling of the chaotic 
part connected to the leads from the entire hierarchical 
part. As the hierarchical region is the semiclassical origin 
of fractal fluctuations, they are not observed. For much 
higher energies only, the hierarchy of partial transport 
barriers would have an outer region with fluxes larger 
than h, leading to fractal conductance fluctuations. The 
inner region of this hierarchy with fluxes smaller than Ti 
has now a smaller phase-space volume. Still, together 
with the regular regions it will cause isolated resonances 
on smaller scales than the fractal fluctuations. Unfortu- 
nately, this energy regime is currently computationally 
inaccessible for the studied system. 

Power laws in the distribution of resonance widths 
and in the variance of conductance increments had been 
observed in Ref. H], apparently reflecting the classical 
dwell time exponent. They were related to the resonances 
below the mean level spacing. For these resonances we 
now have demonstrated that they arise due to hierar- 
chical and regular states. This allows to apply the argu- 
ments of Ref. Jl3| about the resonance width distribution 
of hierarchical states. They lead to the conclusion that 
these apparent power laws come from broad transition 
regions to asymptotic distributions that are unrelated to 
the classical dwell time exponent. 



VI. AVOIDED CROSSINGS 

While for most resonances and corresponding eigen- 
states the parameters rj and T are of the same order of 
magnitude, for some states r) exceeds T by up to 2 or- 
ders of magnitude. This phenomenon can be understood 
as an effect of avoided level crossings in the closed sys- 
tem: In Fig. 0a) we show as an example the dependence 
of the energy of states 5736" and 5737" as a function 
of the parameter M/L for the narrow range 0.10999 < 
M/L < 0.11001, displaying the typical features of an 
avoided crossing. A comparison of the associated Husimi 
representations shows that the states 5736" and 5737~ 
do indeed exchange their character from chaotic to reg- 
ular and from regular to chaotic, respectively, showing 
a superposition at M/L = 0.11. Upon opening the sys- 
tem, the chaotic state couples much more strongly to the 
leads as compared to the regular one. Consequently, in 
the complex energy plane of the scattering system there 
is no longer an avoided crossing. The regular state leads 
to an isolated resonance with an almost linear energy de- 
pendence on M/L and the phase-space signature of the 
regular state (middle row in Fig. |?]b). It closely follows 
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FIG. 7. (a) Energies of states 5736" and 5737" (solid 
lines) showing an avoided crossing under variation of M/L. 
The energy of the only isolated resonance in this energy 
range (dots connected by a dashed line) follows the reg- 
ular state of the closed system. The slight offset in the 
resonance energy is within the systematic numerical er- 
ror of the numerical method, (b) The Husimi representa- 
tions for state 5737" (top row), the scattering state (mid- 
dle row) and for eigenstate 5736" (bottom row) are shown 
for M/L = 0.10999, 0.11, 0.11001 (left to right). For the 
eigenstates one clearly sees the typical exchange of the struc- 
ture while passing the avoided crossing whereas the scattering 
state is not affected. 



the expected energy dependence of the regular state in 
the closed system if it had not made an avoided crossing 
with the chaotic state (Fig. ^|a). 

Another example of an avoided crossing is given by 
state 5801" with an eigenvalue about 0.5A less than the 
resonance position (see lower right corner of Fig. [s]) and 
the state 5802", with an eigenenergy about 1.4 A above 
the resonance energy. Both states show similar Husimi 
representations and have 77 values exceeding T by about 
a factor of ten. 

In all the cases when 77 drastically exceeds V a closer 
look at these eigenstates reveals that they are superposi- 



tions of regular or hierarchical states with chaotic states. 
They are due to avoided crossings and the chaotic part 
leads to a comparatively large value of rj. In contrast, in 
the open system no avoided crossing occurs in the com- 
plex energy plane and the resonance width T is unaf- 
fected. 



VII. CONCLUSION AND OUTLOOK 

We demonstrate a clear correspondence of the isolated 
resonances observed in the transport properties of the 
open cosine billiard to hierarchical and regular eigen- 
states of the closed billiard. We can identify all reso- 
nances with widths less than half of the mean level spac- 
ing. The classification of resonances into a hierarchical 
or regular origin yields numbers in agreement with the 
relative phase-space volumes. On a quantitative level, 
we find a roughly linear relation between the widths of 
the isolated resonances and the weights of the associated 
eigenstates at the part of the boundary where the leads 
are attached. States with unusually large weights can be 
attributed to avoided crossings with chaotic eigenstates. 

We find that the island hierarchy is separated from 
the chaotic part of phase space by partial transport bar- 
riers with fluxes of the order of h. This supports the 
notion that the absence of fractal conductance fluctua- 
tions in the currently accessible energy range is due to 
the quantum dynamical decoupling of the hierarchical 
part of phase space from the chaotic part connected to 
the external leads. 

The simultaneous appearance of isolated resonances 
and fractal fluctuations, beyond the quantum graph 
model studied in Ref. Q, remains to be demonstrated 
numerically or experimentally for a system with a mixed 
phase space. Numerically, the challenge is the observa- 
tion of fractal fluctuations of the conductance, which go 
beyond one order of magnitude p5[ . This requires cal- 
culations with a drastically increased number of modes, 
the use of improved techniques like the modular recursive 
Green's function method P6| , and the search for suitable 
billiard systems where the turnstile fluxes across partial 
barriers are particularly large. Isolated resonances will 
easily appear as soon as the parameter is varied on a suf- 
ficiently small scale. We note, that fluctuations of the 
quantum staying probability, which can be fractal ||, 
cannot show isolated resonances. Similarly, we expect 
no appearance of isolated resonances within the fractal 
fluctuations observed in recent studies |27j, as they are 
unrelated to a classical mixed phase space. 

On the experimental side, fractal conductance fluctu- 
ations have been observed HQ and also isolated reso- 
nances coming from regular regions have been recently re- 
ported [Q . The simultaneous appearance of both types 
including isolated resonances from hierarchical regions re- 
quires to go far enough into the semiclassical regime, i.e. 
to quantum dots with dimensions bigger than 1 /im, as in 
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Ref. M . At the same time the phase coherence time must 
be large enough to resolve isolated resonances of a given 
width and, of course, the parameter, typically a magnetic 
field, must be varied on a sufficiently fine scale. Given 
the experimental limitations it would be helpful if an op- 
timal form for such a quantum dot could be provided by 
theoretical considerations. This seems to be quite diffi- 
cult at present, since the difference in the lithographic 
shape and the actual potential experienced by electrons 
has dramatic consequences on the electron dynamics and 
thus on the scales over which fractal fluctuations and iso- 
lated resonances appear. 
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APPENDIX A: HYBRID REPRESENTATION 
AND RECURSIVE GREEN FUNCTION 
METHOD 

In this appendix we discuss the numerical method to 
determine the scattering matrix S and the time delay r. 
The 5-matrix of a symmetric scattering system can be 
expressed in terms of the Green function G 



(Al) 

(A2) 
(A3) 
(A4) 
(A5) 
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2 / hr , , 
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(A6) 



is the velocity of mode a and 
G af3 (x,x') = f dy f dy'r a (y^)My'^')G + (r,r';E) 

J X J X 1 

(A7) 

is the projection of the retarded Green function 
G + (r, r'; E) onto the local transverse modes 



2 . / any 
■ sin 



W(x) \W(x) 



(A8) 



The Green function can be calculated recursively. Ex- 
panding the Hamiltonian 



H 



d 2 



d 2 



2m \dx 2 dy 2 



(A9) 



in terms of the local transverse modes (A8) and discretiz- 
ing in x direction with a lattice constant a, x = ma, we 
obtain the Hamiltonian in hybrid representation [p9[ 



H h =J2\^rn) (C + 2£ t )(a,m| 

^ {^m, m +\\ a , m ){P, m + 1 \ 

+C+l,ml a > rn + 1 )(A m l) » 



(A10) 



with 



.a/3 



*W(ma) 



tm,m+i =Et C(2/! TO a)0/3(2/; ( m + l)a)dy, 



(All) 



and E t = h 2 /(2ma 2 ). In order to recursively calcu- 
late the Green function associated with , we split the 
Hamiltonian Hff +1 of a system with m = 1, . . . , M + 1 
into two parts, 



H™ +1 = H Q + U, 



Ho = H h 



o 

M 



(A12) 



\a,M + l) (e% +1 + 2E t ) (a,M + l\, (A13) 

f/ = -E(^M+iK M )^ M + 1 l 

a, /3 

+t a J +hM \a,M+l)([3,M\). (A14) 

Dyson's equation, 

G M+1 = G + G Q UG M+1 , (A15) 

then allows to calculate the Green function G M+1 of 
Hjf +1 from G M and 

G Q = (E-H y 1 

= G M +J2\<*,M+ l) 5o M+1 (a, M + 1| , (A16) 

a 

g^ = (E-e« M+1 -2E t )- 1 . 

We start the recursion with M = 1 at the left edge of 
the closed billiard and iterate to the right edge at M = 
Nl = L/a. In order to attach the leads , we again split 
the Hamiltonian according to eq. ( Al2j ), but this time 
Hq contains the Hamiltonian of the closed billiard and 
the semi-infinite leads. In the leads, ef n = (aa/W) 2 and 
Cfro+i = ^apE t - U is the coupling between the billiard 
and the semi-infinite leads. The Green function for the 
leads is known analytically ]30[. 
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The recursion scheme is exact for an infinite number 
Nl of slices and an infinite number Nc of transverse 
modes. For numerical calculations, both numbers have 
to be kept finite. We find that the deviations from the 
asymptotic values for, e.g., the width T of a resonance 
scale as 



T{N Cl N L ) = T + bN c 



(A17) 



with positive numerical coefficients b and c. For our 
choice of parameters, N = 45 transmitting modes in the 
leads, the values Nc — 108 and Nl — 12000 give an 
accuracy of about 1% for the resonance width. The cor- 
rections to the position of the resonance have the same 
functional form as in eq. (A17), however, they can be ei- 
ther positive or negative, depending on the values of Nc 
and Nl . This explains the slight offset of the resonance 
energies with respect to the eigenenergies of the closed 
system seen in Fig. ^a). 
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